Comment on "Characterization of the tunneling conductance across DNA bases" 
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In a recent article, Zikic et al. [Phys. Rev. E 74, 011919 (2006)] present first-principles calculations of the 
DNA nucleotides' electrical conductance. They report qualitative and quantitative differences with previous 
work, in particular with that of Zwolak and Di Ventra [Nano Lett. 5, 421 (2005)] and Lagerqvist et al. [Nano 
Lett. 6, 779 (2006)]. In this comment we address the alleged discrepancies, showing that they come from 
a misrepresentation of our research. Further, we discuss in more detail the issue of geometric fluctuations 
previously investigated by us and raised again in the work of Zikic et al. In addition, we point out erroneous 
comments made by Zikic et al. regarding the use of DFT calculations in transport. 



Recently, Zikic et al. JH report the conductance of pas- 
sivated DNA nucleotides located in between nanoscale gold 
electrodes using density-functional theory (DFT) within the 
known exchange-correlation (xc) functionals. In several 
places throughout their article, they compare their findings 
with previously published results by us J3, Si and conclude 
that there are both qualitative and quantitative differences with 
their work. We point out that Zikic et al. misrepresent the ex- 
isting literature by leaving out important details and, further, 
make comparisons which are at odds with their own approach 
and conclusions. Also, some comments in their work raise 
general questions about the adequacy of static approaches to 
transport and the differences between such approaches. We 
address these issues below. 

Zikic et al. correctly state that the electronic signature of 
nucleotides is strongly dependent on what they call "geometri- 
cal factors". Their work is an explicit demonstration of a well 
known concept: the tunneling current depends exponentially 
on the width of the tunneling barrier, which is here formed 
by the reduced coupling between the electrodes and the nu- 
cleotide. For the case of DNA between electrodes, this means 
that changes in nucleotide orientation modifies their coupling 
to the electrodes, and therefore can drastically change the 
electrical conductance. In addition, if one fixes the sugar- 
phosphate backbone position, the different sizes and geome- 
tries of the bases will cause them to be more or less close to 
the electrodes and therefore cause a difference in their relative 
conductance. Zikic et al. seem to indicate that these conclu- 
sions are qualitatively and quantitatively different from ours. 
Instead, we understood this fact and we stated explicitly in 
Ref. that "[how well the HOMO and LUMO] states couple 
to both electrodes determines the overall magnitude of the rel- 
ative currents [between the bases]". In addition, well prior to 
their work, recognizing the importance of geometrical factors, 
we explored this issue in much more detail by investigating 
realistic structural fluctuations in Ref. Q]. 

Due to the importance of geometry in the conductance (and 
relative conductance) of the different nucleotides, one can 
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not quantitatively compare the results of Ref. 1 with those of 
Refs. 0, In addition, Zikic et al. state that their DFT cal- 
culations give what they call mutually consistent results with 
different exchange-correlation functionals. By mutually con- 
sistent they mean that the ordering of the current magnitudes 
are the same regardless of the xc-functional used. Consistency 
may hold true for the current averaged over all their configu- 
rations and at small bias, but it is obvious from their own work 
that it does not hold true otherwise. For instance, by examin- 
ing either the conductance or the current in their Figures 8, 9, 
10, 11, or 12, one can extract essentially any desired ordering 
in the nucleotides' conductance. 

Thus, we point out that it is incorrect for them to claim 
that our^ results in Ref. |2j are not consistent with our results 
in Ref. 3 based on the change in conductance ordering. This 
claim leaves out crucial facts which are clearly written in our 
papers: from one paper to the other we did change i) the bias, 
ii) the electrode spacing, and Hi) the nucleotide configurations. 
Anyone of these changes can modify the values of the conduc- 
tance, and even the relative conductance. 

Zikic et al. also fail to mention that in our second work, 
Ref.0, we are sampling over more than a thousand nucleotide 
configurations. Thus, one expects that sampling over a non- 
random subset of configurations - as they do in their work 
by only varying one angle - one can obtain different order- 
ings of the conductance. In fact, the alleged discrepancies 
instead highlight and reinforce one of the main conclusions 
of our work (see Ref. 3): in order to successfully sequence 
DNA via transverse electronic transport, each device has to 
be first calibrated by reading a known strand such that the 
current distributions for all four nucleotides can be obtained. 
These distributions are unique to each and every device and 
are determined by the microscopic geometry of the pore and 
electrodes. 

We now want to turn to two important questions: 1) why is 
it that geometric factors are important in the conductance of 
nucleotides, and 2) how can "distinguishability survive aver- 
aging over possible conformations of ss-DNA" (using similar 
words to the question Zikic et al. raise in their conclusions)? 

The answer to the first question is something not stated by 
either Zikic et al. or us. There are two factors that enable one 
to focus mainly on the geometry of the nucleotide-electrode 
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configuration: 1) the HOMO and LUMO energies of the dif- 
ferent bases are close in energy compared to their distance 
from the gold Fermi level. 2) For all four bases, the HOMO 
and LUMO states are delocalized around the base, and thus 
one can substitute the atomic structure of the bases as an ap- 
proximate representation of the spatial extension of the molec- 
ular states. Contrary to the conclusions of Zikic et al., one 
can not say from their results that HOMO and LUMO states 
are less important than geometry, only that when comparing 
molecules of similar factors (1) and (2), that the geometry 
would be the most dominant factor. Further, this leads to a 
very important conclusion: in the case of sensors to detect the 
DNA bases using electrical currents, the nucleotide configu- 
rations have to be at least partially controlled. In terms of a 
nanopore-based device, one way to do this is to use a trans- 
verse electric field induced by the transverse electrodes or by 
an external capacitor across the whole device. [3] Of course, 
there will be other important factors to consider besides ge- 
ometric fluctuations of the nucleotides themselves, including 
the effects of ions. 

The answer to the question on how the distinguishabil- 
ity survives by averaging over possible conformations of 
ss-DNA, can be found by examining, in the context of a 
nanopore-based device, the geometric configurations of DNA 
as it translocates through a pore in the absence of any con- 
trol. This is the basis of our work in Ref. 3. The geometrical 
fluctuations cause the different nucleotides to have large fluc- 
tuations in the value of their current, as shown in Figure Q] 
With no stabilizing transverse field added, the current distri- 
butions for a given set of initial conditions of the different 
nucleotides (calculated as reported in Ref. 3, for an electrode 
bias of 0. 1 V and an electrode spacing of 15 A) span several 
orders of magnitude, and have significant overlap, as shown 
in the figure. These large fluctuations will cause the bases to 
be essentially indistinguishable. However, in the presence of 
a transverse field that is much larger than the driving field, 
the nucleotides can be stabilized, and thus distinguished by a 
relatively modest ensemble measurement 13J. 

We conclude by discussing the differences in using a tight- 
binding (TB) approach compared to DFT calculations for the 
problem at hand. Similar to the DFT approach within any 
available xc-functional, a TB approach has its own limita- 
tions. Nonetheless, in the present context it has a clear ad- 
vantage. In particular, it satisfies two conditions required by 
any method to investigate the relative conductance of the nu- 
cleotides. First, since the coupling (i.e., geometry) is the large 
determining factor in the relative conductance, one needs to 
adequately reproduce the spatial distribution of the molecular 
wavefunctions. Second, the energies of the molecular states 
need to be calculated fairly accurately. For our chosen TB 
parameters, both of these quantities compare well with DFT 
calculations for isolated nucleotides. This, together with the 
reduced computational complexity of TB, allows us to look at 
many different geometric configurations to more realistically 
capture structural fluctuations that would be present in an ex- 
periment. 
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FIG. 1: Probability distributions of currents, for unstabilized 
poly(dX) 15 as the strand propagates through a pore with embed- 
ded electrodes. X is Adenine/Thymine/Cytosine/Guanine for the 
black/blue/red/green curves, respectively. The thin lines show the 
actual current intervals used for the count, while the thick lines are 
an interpolation. 

In addition, there is no reason to believe that DFT, within 
the xc-functionals used by Zikic et al., can represent more ac- 
curately the nucleotide-electrode coupling compared to TB in 
this particular geometry where the nucleotides are not cova- 
lently bonded to either electrodes. [4] Indeed, the exchange- 
correlation functionals employed in Zikic et al.'s work do 
not include the long-range van der Waals interactions that are 
present in this weak coupling case. The fact that the two dif- 
ferent xc-functionals employed by Zikic et al. show order of 
magnitude differences in the conductance may be a result of 
this problem. 

Finally, Zikic et al. state that 'As far as the self-consistency 
of the electron transport is considered, this leads to a pro- 
cedure that is asymptotically exact in limit of zero electric 
bias". This is a misconception about static approaches to 
transport, and is not correct even if one had the exact static 
xc-functional. Two of the present authors (MZ and MD) have 
shown [5] that even in the limit of zero bias, with the inclusion 
of self-consistency, the current obtained using static DFT does 
not include dynamical many-body effects which can only be 
captured by using time-dependent approaches such as time- 
dependent DFT. [6] In Ref. 5 we have evaluated these dynam- 
ical corrections specifically for the local density approxima- 
tion functional, however, the statement is true regardless of the 
static xc-functional chosen: no static xc-functional (even the 
exact one) can fully capture the true non-equilibrium nature of 
transport problems. Incidentally, these dynamical many-body 
effects are also absent in the TB static approach to transport 
we have employed. 
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